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Abstract 

Complementary Neutron Spin Echo and X-ray experiments and Molecular Dynamics simulations have 
been performed on difluorotetrachloroethane (CFC12-CFC12) glassy crystal. Static, single-molecule re- 
orientational dynamics and collective dynamics properties are investigated. The orientational disorder is 
characterized at different temperatures and a change in nature of rotational dynamics is observed. We show 
that dynamics can be described by some scaling predictions of the Mode Coupling Theory (MCT) and a 
critical temperature T c is determined. Our results also confirm the strong analogy between molecular liquids 
and plastic crystals for which a-relaxation times and non-ergodicity parameters are controlled by the non 
trivial static correlations as predicted by MCT. 
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I. INTRODUCTION 



Understanding of the mechanisms leading to the glass transition has been a subject of intense 
research in recent years and it is still a matter of numerous theoretical and experimental studies Q]. 
A fundamental question particularly concerns the microscopic description of the precusor coop- 
erative mechanisms which develop in the high temperatures pico-nanosecond regime (ps-ns) far 
above the glass transition temperature T g . It is now commonly suspected, mainly from theoretical 
and numerical results, that dynamical behaviours of glass-forming systems could be described by 
involving two characteristic temperatures [2]: i) (Ts or T x are also seen) corresponding to the 
onset of slow dynamics where non-exponential relaxation and non-Arrhenius behavior emerge 
First appearance of rotational heterogeneity are also reported O]. The role of T4 in the vitrica- 
tion process is still not completely understood but it is often considered as the temperature below 
which dynamics are influenced by the potential energy landscape 0]. ii) T c , the critical tem- 
perature predicted by the Mode Coupling Theory MCT [4]. In the last two decades, enormous 
work has been devoted to this theory and it is so far the only one which provides a microscopic 
description of simple liquids. The intrinsic basis of MCT states that the behaviour of any time- 
dependent correlators describing the dynamics of the system is only controlled by its static density 
correlator i.e S(Q) and predicts the existence of a critical temperature T c indicating a bifurcation 
from ergodic to non-ergodic dynamics. The role of S(Q) and the scaling properties predicted by 
MCT approaching T c have been successfully validated in numerous experiments and molecular 
dynamics (MD) simulations [5]. However, no sharp transition is observed at T c itself owing some 
activated processes not included in the theory which restores ergodicity. Thus, T c is also associated 
to a dynamical crossover to a 'landscape-dominated' regime in which diffusion processes can be 
described in terms of thermally activated jumps. 
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Most of the glassformers are molecular liquids and accordingly possess both translational and 
rotational degrees of freedom (TDOF and ODOF). Approaching the glass transtion, the coupling 
between TDOF and ODOF is far from being understood and the precise role of the ODOF remains 
unclear as demonstrated by the so-called low-temperature translation-rotation paradox ||6|, |2U. 
MCT has been initially developped to provide a microscopic understanding of simple atomic 
liquids. However, general predictions have also been successfully validated for molecular com- 
pounds. Recently, some extensions of MCT called molecular mode coupling theory (MMCT) o] 
have been proposed to take ODOF into account: one diatomic probe molecule in an atomic liq- 

i n n 

uid [[8J], liquids made of linear molecules |9] or water H lQfl to cite only a few. The authors have 



particularly shown that some of the basic predictions of MCT still hold owing TDOF/ODOF cou 



1211 performed on orthoterphenyl (OTP), 



pling. As revealed by recent MD investigations 
coupling of the rotational dynamics to the center-of-mass motion can be complex. 

An alternative to the glassforming molecular liquids is offered by at least two partially disor- 
dered systems which provide an interesting solution in order to investigate only ODOF during the 
freezing process : 

i) The so-called orientational glasses such as Ar!_ x (N 2 ) x or (KCN) x (KBr)i_ x mixed crys- 
tal where rotational motions are frustated by a quenched disorder externally imposed by a dilu- 



tion 1 13] K. H. Michel and coworkers B14I1 developped a MCT for these systems which involved 
the interplay or rotational-translational and rotational- strains coupling where strains result from 
the substitutional disorder. 



18], orthocarborane [19] 



ii) Glassy crystals |15] such as cyanoadamantane |16], ethanol |17 

I I 

or cyclooctanol [20] which are true rotational analogs of canonical liquid glassformers. These 
systems are generally composed of globular molecules and exhibit a partially disordered phase, 



called plastic, which possesses a crystalline translational order and a dynamical orientational dis- 



order. They can be rapidly cooled from the plastic phase and exhibit a step of the specific heat at 
the calorimetric temperature T g or a non-Arrhenius behavior of the rotational relaxati on time s. 



22, 



23], 



No microscopic theory has been developped for plastic crystals so far. However, in (21 
we have particularly shown from NMR and Raman experiments, and MD computer simulations 
that some predictions of the idealized version of the MCT (critical temperature T c and time scal- 
ing laws) were able to describe relatively well rotational dynamics of different plastic crystals. It 
should be noted that the other remarkable temperature Ta was also identified. These intriguing pre- 
vious results call for new investigations to clarify the similarity between slow dynamics behaviour 
of plastic crystals and molecular liquids and in particular to discuss the role of the static correla- 
tions S(Q). In this paper, we present for the first time results of Neutron spin-echo (NSE) and 
X-ray experiments and MD simulations of the glassy crystal difluorotetrachloroethane (DFTCE). 



This compound is composed of simple molecules CFC1 2 — CFC1 2 close to dumbbe 
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Is extensively 



26]. DFTCE has 



used in MD calculations as prototype of molecular glass-formers liquids [24 
been experimentally widely studied and presents a rich variety of unusual properties. DFTCE ex- 
hibits a glass transition of the overall rotation of the molecules at T g = 86 K H27I1 . Two additional 
heat anomalies are found at 60 and 130 K and associated respectively to a sub-T 9 (3 process and the 
freezing of the transformation between trans and gauche conformation of the molecule. Changes 



concerning the nature of dynamics in this system has been reported from NMR experiments [29], 



Brillouin and dielectric spectroscopy 112 811 . In this latter investigation, a crossover from individual 
to correlated rotational motions was particularly suggested. 

II. EXPERIMENTS AND DETAILS OF THE SIMULATION 



The samples used in both X-ray and NSE experiments consist of commercial 1,2- 
difluorotetrachloroethane purchased from Sigma- Aldrich and purified by melting zone. 
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The X-ray diffraction experiment was performed on a laboratory diffractometer equipped with 
an INEL curved sensitive detector CPS 120 made up of 4096 channels and allowing us to record si- 
multaneously X-ray diffraction on a 29 range of about 120 deg. The incident beam was monochro- 
matized with a bent of quartz crystal which selects the K ai wavelength of a Cu X-ray tube 
(A = 1.54056 A). Powder DFTCE was introduced in Lindemann glass capillary of 0.5 mm in 
diameter. It was mounted on the axis of the diffractometer and rotated during the experiment in or- 
der to reduce the effect of preferential orientation. The correspondence between the channel num- 
ber and the 29 angle was determined by a preliminary X-ray experiment of cubic Na 2 Ca 2 Al2F 14 



(NAC) 113011 with a cubic spline interpolation between the Bragg peaks of NAC. Low temperatures 
were achieved by placing the sample into a gaseous nitrogen flow at the required temperature. 
Using this procedure, the sample is quenched within a few seconds, and its temperature can be 
controlled within 0.5 K between 1 10 K and 300 K. In the present experiment, the X-ray pattern 
was recorded at 140 K during 3600 s for a sample quenched at 2K/min. 

The NSE experiments were performed on the INI 1 spectrometer at the Institut Laue-Langevin 
(ILL), Grenoble, France. We used the multidetector version of INI 1 which allows simultaneous 
measurements at different scattering angles. A setup with an incident wavelength of A = 5.5 A 
(AA/A ~ 16 %) gave access to wave numbers approximately between 1.28 and 1.72 A -1 and 
we covered the time range 8 ps - 1 ns. We made measurements of the coherent signal of poly- 
crystalline systems filled into a flat Al container from T = 130 to 220 K in steps of 10 K. For 
normalization, a resolution scan of the elastic scattering was realized at T = 5 K in the glassy 
phase and all data were divided by the polarization at this temperature. Two samples of different 
mass identified as A (m ~ 0.33 g) and B (m ~ 2 g) were investigated. Although preliminary 
measurements showed that transmission of sample A was better, 93% to be compared to 77% 
obtained for sample B, this latter was prefered for most of the measurements since its scattering 



signal was more intense. Unfortunately, only data obtained from sample A are available at T = 
130 and 140 K. 

MD calculations were performed on a system of iV = 686 (7 x 7 x 7 bcc crystalline cells) 
molecules. Each DFTCE molecule CFC1 2 — CFC1 2 is described by its 8 atoms and considered as 
a rigid unit. They interact through a Buckingham short range atom-atom potential (see parameters 
in table 0]) and the electrostatic contributions have been neglected since the DFTCE molecule 
possesses a weak dipolar moment. No structural or dynamical change has been found for a system 
where electrostatic interactions are taken into account. It should be noted that both structure 
and dynamics were found in good agreement with experimental results (see figures [T] and |3]>. 
Newton's equations of motion were solved with a time step of At = 5 fs. We worked in the NPT 
statistical ensemble with periodic boundaries conditions where the simulation box is allowed to 
change in size and shape. MD simulations were done for a sample corresponding to the DFTCE 
orientationally disordered phase at 14 different temperatures from T = 130 to 260 K in steps of 
10 K. It should be mentioned that our very simple DFTCE model allows us to perform very long 
MD runs of about 50 ns. 

In this paper, single-molecule reorientational dynamics of DFTCE will be obtained from the 



first two self angular correlation functions C/ = i 2 Il24ll defined as: 

C ^) = Jj E (Pi{ui{t)-uM)) (1) 

i=l,N 

where Pi is the /-order Legendre polynomial and U{ is a unit vector directed along the C3-symmetry 
carbon-carbon molecular axis of the DFTCE molecule i (see Fig.[l]). From C/ = i j2 , we can obtain 
clarification on the nature of dynamical changes and anisotropy of the molecular motions. Ex- 
perimentally, Ci(t) can be directly measured in dielectric relaxation and C 2 (i) in Raman scatter- 
ing Q31 L Cj = i j2 can also be related to the informations obtained from NMR relaxation measure- 
ments ii2 ih . 



In order to investigate collective dynamics, the intermediate scattering function as it can be 
classically obtained from coherent neutron scattering experiments will be mainly considered: 

S(Q,t) = (p (t)p (Q)) (2) 

where pMt) is the time-dependent density correlator: 

Pq(^) = ^2b a exp[iQ.f a (t)} (3) 

a 

where the sum is over all the atoms a of the system. b a and f a are the coherent scattering length 
and the position of the a atom respectively. An average over isotropically distributed Q-vectors 
having the same modulus Q is performed in order to obtain S(Q,t) for a poly crystalline sample. 



In general, S(Q,t) can be expressed as IB2I1 



S(Q,t) = S c (Q) + S d (Q,t) (4) 
where S C (Q) = \(pq)\ 2 is the coherent elastic scattering and 

S d (Q,t) = (SpfflSpdO)) (5) 

where 5pq(t) = Pq*(t) — (p^) is the fluctuation of the time-dependent density operator. For plastic 
crystals, S(Q,t) is identical to S d (Q,t) except for Q-vectors corresponding to Bragg peaks for 
which the long time limit of S(Q, t) reaches the non-zero value S C (Q). This behaviour is found 
at all temperatures and it is associated to the crystalline order of the TDOF and not the freezing 
process of the ODOF. Similarly to the intermediate scattering function, the static correlator can be 
decomposed into two components S(Q) = S C (Q) + Sd(Q) where Sd(Q) is the diffuse scattering: 

S d (Q) = (\Spq\ 2 ) = (\Pq\ 2 }-\(Pq}\ 2 (6) 
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HI. RESULTS AND DISCUSSION 



A. Static 

Owing the high rotational disorder, plastics crystals yield only a few diffraction peaks but ex- 
hibit intense and highly structured diffusive scattering which provides useful information to char- 
acterize rotational motions and most probable orientations IB2I1 . Fig. [l] shows both numerical and 
experimental static structure factors obtained in the Q-range [0.3 — 6.5] A -1 . The diffraction pat- 
tern of the plastic phase presents eight Bragg peaks localized at 1.272, 1.799, 2.203, 2.543, 2.844, 
3.114, 3.596 and 3.816 A" 1 which are well indexed as a bcc structure having a cell parameter of 
6.985 A -1 . The cell parameter is found 2% smaller from MD simulation. This result confirms 
the preliminary X-ray diffraction investigation of this crystal performed by Kishimoto et Al [27]. 
The elastic and the diffuse scattering have been extracted from the MD study and are found in 
good agreement with the X-ray experiment: maxima of S C (Q) coincide with the position of the 
experimental Bragg peaks and maxima of Sd(Q) with the broad bumps of the diffuse scattering 
at Q ~ 1.3, 2.34, 3.36 and 4.8 A -1 . It should be mentionned that an additional diffuse scattering 
bump is observed experimentally close to Q ~ 1.6 A -1 . This feature which is also present in the 
pattern of the liquid phase (see Fig.[TJ mainly derives from the scattering contribtuon of the glass 
capillary sample holder. Moreover, a strong similartiy of the X-ray pattern is observed between 



the liquid and plastic phase. Such remarkable feature has been already reported for ethanol [33] 

n 

or CBr 4 [32]. The first sharp diffraction peak of the liquid particularly corresponds to the first 
and most intense Bragg peak of the plastic phase. Broad bumps of the diffuse scattering are found 
almost identical in both phases and merge at the largest Q-vectors since the same intramolecular 
structure is probed. This result could suggest a strong analogy between the supercooled liquid and 
the supercooled plastic crystal of the same compound in the framework of the MCT which predict 



that dynamics are controlled by non trivial static correlations. This question will be adressed in an 



other paper Il34ll 



The orientational distribution of the DFTCE molecules can be described with the probability 
distribution function f(u) where u = (x, y, z) is a unit vector along the carbon-carbon molecular 
axis of DFTCE (see Fig.[l]). Since the site symmetry of the crytal is O h , f(u) may be expanded as 
a series of cubic harmonics Ki 113 511 : 



\ i=4,6,8,... / 



where K 4 = (21/16) 1/2 (5Q - 3), K 6 = (13/128) 1/2 (462 J R + 21Q - 17) and K 8 = 
(561/1024) 1 / 2 (65g 2 - 208i? - 94Q + 33) with Q = x 4 + y 4 + and R = x 2 y 2 z 2 . The (Ki) 
coefficients are the averages of the corresponding Ki. K 4 , K 6 and K$ have been calculated from 
MD simulations at different temperatures and are reported in table HH This table clearly shows 
no strong ordering of the DFTCE molecules along the prefered [1 1 1], [1 10] or [100] directions of 
the cubic lattice consistent with the high rotational disorder of plastic phases. However, a slight 
preference for [111] directions is observed. Upon lowering the temperature, we also found an 
increase of the orientational probability along [111] unlikely to [110] or [100]. These results have 
been also confirmed from different snapshots of individual molecule trajectories as displayed in 
FigEl 



B. Single-molecule reorientational dynamics 

Single molecule reorientational dynamics have been investigated from the self angular corre- 
lation functions Ci=i t 2 defined in Eq. d Figure El shows (^(t) time correlation functions for all 
investigated temperatures from 130 K to 260 K. Ci(t) time correlation functions, not displayed in 
the present paper, shows a very similar behaviour to the C 2 (t) functions. At high temperatures, 

9 



C 2 (£) rapidly decays zero and can be described by a simple exponential shape. When lower- 
ing the temperature, dynamics slow down and a two-step process is seen as already observed in 



chloroadamantane plastic crystals 1121 . 



23H and in supercooled liquids B24I1 . This process is often 
interpreted in term of a— (5 relaxations which will be carefully analyzed in the next section for col- 
lective dynamics in the framework of MCT. Reorientation times t\ and t 2 have been extracted from 
the time it takes for their respective self angular correlation functions C/ = i j2 to decay e _1 of their 
initial values. The ratio ti/t 2 is displayed in Fig. El For uncorrelated small angular steps motion, it 
can be shown that Ci (t) functions follow an exponential decay exp[— t / n] where 1/tj = 1(1 + l)D r 
and D r is the rotational diffusion coefficient. At high temperatures, ti/t 2 is close to the value of 3 
in good agreement with this simple diffusion model. Upon lowering the temperature, ti/t 2 clearly 
deviates from 3 and approaches 1 . In [36], we have shown that this behaviour indicates a crossover 
to large angle tumbling motions. Similar features have been recently shown from MD simulations 



of supercooled liquid CS 2 113711 . This result is also confirmed from snapshots of single-molecule 
trajectories in Fig.|2| Thus, low temperature rotational dynamics can be seen as tumblings between 
prefered orientations along [111] directions of the cubic cell. This result also proves that simple 
exponential decay for Ci(t) function is not granted at low temperature and that non-exponential 
relaxations occur. 



C. Collective dynamics 

The coherent intermediate scattering function S(Q, t) has been calculated from MD at different 
temperatures and different Q wave- vectors and it is compared with the NSE data in Fig. |3] at the 
arbitrary wave-vector Q = 1.54 A -1 . A good agreement between experimental and numerical 
S(Q, t) is found. However, some discrepancies are particularly observed at high temperatures at 
the longest Fourier times for which the NSE S(Q, t) does not decay completely to zero. Since no 
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elastic contribution is expected outside the Bragg Peaks, this feature is most likely due to uncer- 
tainties in the background substraction. At low temperatures, at T = 130 and 140 K, difference 
between experiment and simulation are also seen. Those are also the data which show the largest 
error bars. One explanation could be the use of a different sample (see experimental details) pro- 
viding a less intense scattering signal. 



Upon cooling, S(Q,t) displays a decay classically observed in MD simulations [121 and in 
NSE experiments of glass-forming molecular liquids such as glycerol [38] or CKN 10911 . For 
DFTCE glassy crystal, the long time relaxation corresponding to the overall rotations of molecules 
is separated from short time regime by a plateau-like region. This latter can be associated to the 
orientational trapping of molecules, the analogue of the cage effect in liquids [40], as already 
demonstrated in [36] for others plastic crystals. In the long time relaxation regime, we defined the 
characteristic time ru e of the rotational dynamics as the time it takes for S(Q, t) to decay from 
1 to 1/e at the wave- vector Q = 1.54 A -1 . This arbitrary definition is chosen in order to fully 
compare both NSE and MD data. Owing the limited Q and time domain investigated in the present 
NSE experiment, other definitions such as fits with stretched exponential of Kohlrausch or time 
measured to decay from 1 to 0.1 can not be employed. An alternative definition will be particularly 
used in the following section concerning the MCT analysis of the MD data. 

Ti/ e obtained experimentally and numerically are displayed in Fig. |4j We found that both 



NSE and MD characteristic times T\u are consistent with previous NMR 12911 and dielectric 
investigations and allows us to extend DFTCE data at high temperatures over 3-4 decades. It 
clearly confirms the non-Arrhenius behavior of the rotational dynamics. For comparison, we also 
plotted in Fig. 13 the Vogel-Fulcher-Tammann (VFT) law tvft = Too gxp[1/Kvft-(T/Tvft — 1)] 
where Tvft is the temperature of apparent divergence of tvft and Kvft measures the kinetic 
fragility. Parameters K VFT = 0.13 and T VFT = 70 K were determined from dielectric exper- 
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iments of plastic DFTCE B28H in the temperature range [95-120] K. From MD simulation, we 
obtained Kvft — 0.20 and Tvft — 86.9 K for t\u in fair agreement with the dielectric experi- 
ment. A perfect agreement can not be expected owing the Q-dependence of the neutron scattering 
technique and the different temperature domain investigated. 



In J23 



3611 - we suggested that some common microscopic mechanisms, relatively well de- 
scribed by MCT, are involved in both orientationally disordered crystals and molecular liquids. 
Therefore, in the following, we analyse our data in the MCT framework. In Fig.|6l S(Q, t) is plot- 
ted against the rescaled time T\i e previously determined. We find that one of the main predictions 
of MCT i.e the time-temperature superposition principle (TTSP) holds relatively well both for 
NSE and MD data. We show that the curves overlap for the large time part and collapse gradually 
onto a master curve upon cooling. This feature is made particularly clear since we have developed 
a realistic MD model of DFTCE which allows us to access data in the plateau region of S(Q, t) 
which was not possible in the present NSE experiment. At T =130 K, S(Q, t/ r i/e) obtained from 
MD deviates from the master curve. This indicates that MCT ceases to be valid at that temperature. 
This result is confirmed from snapshots and trajectories of target molecules, as already shown for 
other plastic crystals in 



36], which reveal that at the lowest temperatures rotational motions 
correspond to activated tumblings of the overall molecules between preferred molecular orienta- 
tions i.e the so-called hopping processes not taken into account in the ideal version of MCT (see 
section IlII All . Moreover, this deviation also allows us to roughly estimate the critical temperature 

T 

± c- 

Above T c , MCT predicts that the late (3 regime or the early a relaxation is described by the 
following power law (going beyond first order): 



t) = %-h$.(t/Ty+h$.(t/T)» (8) 

where f% is the so-called nonergodicity parameter associated to the plateau height and r = t (T — 
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T c ) 7 with 7 = 1 /2a + 1/26. The parameters to>7 and b are temperature and q independent and 
related via 7 = l/2a + 1/26 and T 2 (l - a)/r(l - 2a) = T 2 (l + 6)/r(l + 2b) where is the 
gamma function. The first two terms of Eq.[5]correspond to the classical von Schweidler law and 
the last term is a second order correction. First, we fixed the exponent b by performing individual 
fits of the different correlators S(Q, t) in the late (3 regime using Eq.[8]at T = 140 K expected to 
be the lowest temperature where MCT predictions are valid (see Fig. We used the following 
free parameters: b, /§, & 1 ) = h^(Q).T- b and h {2) = hW{Q).T~ 2b . Best results are obtained 
for b = 0.59 which corresponds to a = 0.31 and 7 = 2.45. Then, using this value of 7, scaling 
properties in the a relaxation regime have been checked using two functions Ri(T) and R2(T) 
defined as: 

T- lh (T)-r lh (Tl) 
Ri(T) = _ i /-, . \ _ lAv/ \ (T2 - Tl) + Tl (9) 



T, 



- X '\T2)-To Xh {Tl) 



Q "J >Q 

R 2 (T) = — 3 — — (T2 - Tl) (10) 

where T\ and T 2 are two arbitrary low temperatures and tq(T) is a characteric time of the a 
relaxation defined as the time it takes for S(Q,t)/S(Q,t = 2 ps) to decay from 1 to 1/e. This 
alternative definition eliminates the short-time dynamics and ensures that the (3 relaxation does 
not pertub the a-relaxation time i.e. tq(T) is not influenced by other parameters f Q , h^(Q) or 
b. Therefore, it avoids a tedious fitting procedure with stretched exponential of Kohlrausch. It 
is also assumed that any characteristic times belonging to the a regime show asymptotically the 
same temperature dependence tq ~ r defined in Eq. [8j For all Q vectors, it can be shown that 
Ri (T) = T and R 2 = T — T c in the temperature range where MCT predictions are valid. In Fig.0 
Ri and R 2 scaling funtions are displayed as function of temperature at different Q wave vectors. 
They clearly indicate a temperature range [140-180] K where all data obey an universal master 
curve. This result is particularly significative owing the wide Q-range investigated and confirms 
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the choice for the exponent b. Both R\ and R 2 allow us to estimate two remarkable temperatures: 
i) Ta — 180 K corresponding to the crossover from simple liquids dynamics to slow dynamics. As 
observed in previous MD of molecular glass-forming liquids W24I1 . polymer melts 114 ill or plastic 
crystals [23], this temperature seems very dependent on the chosen correlator S(Q,t) i.e on the 
space range probed at the wave- vector Q. ii) the critical temperature predicted by MCT T c ~ 128 
K. The superposition of the different correlators seems better for R\ than R 2 . Thus, contrary to 
MCT predictions, this reveals a weak Q-dependence of T c of about 3 K. In order to check the limit 
of this approach, the same work was also performend on the NSE data and results are displayed as 
inset in Fig|6l It should be mentionned that tq is not accessible in our present experimental data 
and we had to use T\u which in intrinsically worse. T\/ e is influenced by the /3-regime. Moreover, 
b obtained from MD was employed. Despite these two drawbacks, reasonable results are obtained 
for both Ri and R 2 for the narrow Q-range investigated in fair agreement with the MD data and 
clearly mmotivate measurement at shorter times using time-of-flight experiment to describe the 
plateau regime. 

A fundamental property of MCT or MMCT stresses that dynamics are driven by the static 



density correlator S(Q). This 



has 



:>een confirmed in several numerical and experimental studies 



for simple W42I1 or molecular 11121 14311 glass-forming liquids. Since a good agreement between 
experiment and simulation was found for the elastic and the diffuse structure factor (see Fig. [T]), 
we checked the Q-dependence of the a-relaxation time tq, the non-ergodicity parameter /q and 
the total prefactors h^(Q) and h {2 \Q) at T = 160 K. Results are shown in Fig. |7] Clearly, 
we found that tq, fq oscillate in phase with Sd(Q) and S C (Q) while the amplitudes h^(Q) and 
h^ 2 \Q) oscillate in phase with Sa{Q) and out of phase with S C (Q). It should be noted that the Q- 
modulation of the non-ergodicity parameter and the total prefactor strongly ressembles to results 
obtained by Franosch et Al. for a hard-sphere colloidal suspension model 114411 . 
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IV. CONCLUSION 



In this paper, we have investigated static and dynamical properties of difluorotetrachloroethane 
glassy crystal using complementary X-ray and NSE experiments and MD simulations performed 
for a realistic molecular model. This study clearly confirms that supercooled molecular liquids and 
supercooled orientationally disordered crystal share many properties in the ps-ns regime where 
some precursor features (non-exponentiality or non-Arrhenius behavior) of the glass transition 
emerge. We particularly found that some MCT equations developped for atomic liquids provide a 
good decription of rotational dynamics: some MCT predictions such as the superposition principle 
or the power law (T — T c )~ 7 dependence (T c ~ 128 K )of the a -relaxation time caracterizing col- 
lective dynamics hold well. We also demonstrated that the a -relaxation time tq, the non-ergodicity 
parameter /q and the amplitudes hi(Q) were modulated by the static correlation consistent with 
the MCT predictions. In addition, the other remarkable temperature T A was detected at about 180 
K for both self and collective dynamics. Therefore, in order to fully validate the MD results, it 
will be of great interest: i) to complete the present NSE study over a wider Q-range and for longer 
Fourier times in order to reach the estimated T c and ii) to make measurement at shorter times using 
time-of-flight experiment to describe the plateau regime. 
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Tables captions 

TABLE I: Simulation coefficients taken from literature for Buckingham potential 4>(r) = Aex.p(—p.r) — 
C/r e . 

TABLE II: Orientational order parameters (K4), (K§) and (K$) at different temperatures. Values for pure 
[111], [110] and [100] ordering are also indicated. 
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TABLE 1 



Site-Site A (kJmol- 1 ) p (A" 1 ) C (kJmol^A 6 ) 



c - 


c 


226307 


0.288 


2420 


C - 
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F 
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1 96747 


260 


1 168 


c- 


CI 


390940 


0.284 


3864 


F- 
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171038 


0.237 


565 


F- 


CI 


320883 


0.258 


1808 


Cl 


-CI 


586389 


0.284 


5798 
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TABLE 2 



Temperature (K4) (K 6 ) (K 8 ) 



140 


-0 647 465 211 


160 


-0 638 409 191 


180 


-0 623 368 173 


200 


-0 596 319 151 


220 


-0 577 291 134 


240 


-0 549 257 1 20 


260 


-0.515 0.223 0.105 


[111] 


-1.528 2.266 0.877 


[110] 


-0.573 -2.071 1.665 


[100] 


2.291 1.275 2.961 



Figures captions 



FIG. 1 : Experimental X-ray powder diffraction pattern in the liquid phase at T = 305 K(top) and the plastic 
phase at T = 140 K (bottom) as function of the wave-vector Q. The calculated coherent elastic scattering 
S C (Q) (dotted line) and the calculated diffuse scattering Sd(Q) (dashed line) in the plastic phase are also 
indicated for comparison. The DFTCE molecule is displayed in inset. 

FIG. 2: Trajectories (top) of the unit vector u(t) = (x(t),y(t), z{t)) of one target molecule over a 1 ns 
run at 140 K (left) and 240 K (right). Horizontal lines indicates prefered orientation of u along cubic [111] 
directions. Ratio (bottom) of reorientational correlation times t% and T2 obtained from MD simulation as 
function of temperature from 130 to 260 K The reorientational correlation function C2(t) as function of 
time at different temperatures from 130 to 260 K is given in inset. 

FIG. 3: Coherent intermediate scattering function S(Q,t) at Q = 1.54 A 1 obtained from Neutron Spin 
Echo experiments (full circle) and MD computer simulations (solid line) at T = 130, 140, 150, 160, 180 
and 200 K. 

FIG. 4: Characteristic time t\u obtained from NSE experiment (full circle) and MD simulation (open circle) 



at Q = 1.54 A 1 as function of the temperature. Data extracted from NMR (square) [29] and dielectric 



experiment using a VFT fit (gray line) [28] are also indicated. 

FIG. 5: S(Q, t) obtained from NSE experiment at Q = 1.54 A -1 plotted against the scaled time i/ T i/e at 
different temperatures from T = 140 to 190 K. S(Q, t) calculated from MD simulation at same Q wave- 
vectors is displayed at different temperatures from 130 to 220 K. 
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FIG. 6: Scaling functions R\ and R 2 at different wave-vectors Q ranging from 0.7 to 3.9 A 1 as function 
of the temperature. R\ and R2 are defined in Eq. [9] and [^respectively using exponent 7 = 2.45 extracted 
from the MCT analysis in the (3 regime and T x = 140 K and T 2 = 170 K. R 1 (T) = T and R 2 (T) = T - T c 
are expected from MCT predictions. R\ and R 2 extracted from NSE experiment using characteristic times 
r 1/e , Ti = 150 K and T 2 = 190 K are given in inset. 

FIG. 7: Q-dependence of different parameters at the temperature T = 160 K. a) The a-relaxation times tq 
extracted from the Sd{Q,t) correlators, b) and c) The non-ergodicity parameter /q and the total prefator 
h^ l \Q) obtained from a fitting procedure of S(Q,t) using a von Schweidler law as derived from MCT 
including second order corrections (see Eq.[8l>. hP>(Q) is not shown but exhibits the same Q-dependence 
as ^'{Q). d) The diffuse scattering Sd(Q)- First Bragg peaks positions clearly identified from MD are 
also indicated with vertical lines. 
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FIGURE 2 
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FIGURE 4 



GO 



10 

io 7 
io e 
io 5 
io 4 



b* 3 

10 



10 



10 



10 



□ NMR 
O MD 
• NSE 

VFT (dielectric) 



□ 



o 



o 

□ 



□ 



□ 



□ 



□ 




4 



6 7 8 

1000/T 



10 



27 




28 



FIGURE 6 
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